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摘要 : 超 导 同 步 回 旋 加 速 器 具有 超 小 型 结构 ,造价 更 低 ， 更 适合 应 用 于 医院 环境 使 用 的 质子 
治疗 系统 。 相 比 于 现在 的 质子 治疗 系统 , 应 用 了 超 导 同 步 回旋 加 速 器 的 质子 治疗 系统 单 次 进 
行 肿 瘤 治 疗 的 费用 更 低 ， 具 有 巨大 的 应 用 前 景 ， 因此 取得 许多 加 速 器 研发 机 构 和 企业 的 高 度 
重视 。 再 生 引出 系统 是 超 导 同 步 回 旋 加 速 器 的 关键 系统 , 也 是 超时 同步 回旋 加 速 中 设计 的 关 
键 和 难点 。 本 文 针对 一 台 质子 治疗 用 230MeV 质子 超 导 同 步 回 旋 加 速 器 开展 了 再 生 引 出 系 
统 设 计 。 由 于 再 生 引出 系统 与 主 磁场 紧密 相关 ,文章 第 一 部 分 先 开 展 了 主 磁铁 的 设计 ， 并 给 
出 相关 的 动力 学 计算 结果 。 现 有 的 回旋 加 速 器 设计 程序 并 没有 针对 同步 回旋 加 速 器 再 生 引 出 
系统 的 设计 功能 , 故 本 文 开发 了 能 够 进行 再 生 引 出 系统 设计 的 程序 。 设计 结果 显示 230 MeV 
质子 超 导 同 步 回 旋 加 速 器 再 生 引出 系统 的 设计 满足 引出 要 求 , 可 以 为 超 导 同 步 回旋 加 速 器 其 
他 系统 的 设计 及 工程 建造 提供 参考 。 

关键 词 : 超 导 同 步 回 旋 加 速 器 ， 超 小 型 结构 ， 再 生 引 出 系统 设计 ， 主 磁铁 设计 ; 


0 引言 : 
传统 的 放疗 手段 利用 光子 束 辐 照 肿瘤 组 织 , 然而 光子 束 辐 照 人 体 时 在 人 体 中 沉积 包 


而 在 肿瘤 区 沉积 的 能 量 更 少 ， 因 此 传统 放疗 手段 存在 严重 副作用 。 质 子 束 在 人 体 中 沉积 

的 能 量 与 速度 的 平方 成 反比 ， 主 要 沉积 在 射程 末端 形成 剂量 “布拉格 峰 ”， 因 此 可 以 控 

制 质子 的 能 量 和 轨迹 ， 使 质子 放射 剂量 精准 沉积 在 肿瘤 区 ， 可 极 大 降低 副作用 。 
同步 加 速 器 、 等 时 性 回旋 加 速 器 和 同步 回旋 加 速 器 是 目前 质子 治疗 系统 采用 的 主要 
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加 速 器 类 型 。 这 三 者 中 ， 同 步 加 速 器 的 尺寸 相对 较 大 、 造 价 高 ， 限 制 了 其 的 推广 
使 用 ; 等 时 性 回旋 加 速 器 由 于 结构 紧凑 ,造价 较 低 ,是 目前 商用 质子 治疗 系统 的 
主流 加 速 器 类 型 ; 同步 回旋 加 速 器 工作 在 脉冲 状态 ,其 平均 流 强 要 远 小 于 回旋 加 

速 器 。 但 是 同步 回旋 加 速 器 磁铁 结构 简单 ， 当 采用 超 导 技 术 时 ， 其 尺寸 和 重量 会 
得 到 大 幅度 减 小 , 使 得 其 可 以 直接 安装 在 旋转 机 架 上 ， 从 而 大 大 降低 束 流传 输 系 
统 的 建造 成 本 , 使 其 成 为 质子 治疗 用 加 速 器 的 有 力 竞争 者 。 国内 外 开展 了 很 多 关 
于 超 导 同 步 回旋 加 速 器 的 设计 和 研制 工作 ， 其 中 Meviong 已 经 实现 了 超 导 同 
步 回 旋 加 速 器 质子 治疗 系统 的 商业 化 生产 。 

束 流 引 出 系统 是 超 导 同 步 回 旋 加 速 器 的 关键 系统 , 也 是 超 导 同 步 回旋 加 速 器 
设计 的 关键 技术 难点 。 一 方面 ， 超 导 同 步 回旋 加 速 器 加 速 粒 子 的 过 程 中 ， 射 频 电 
场 频 率 需 要 根据 粒子 回旋 频率 的 改变 而 改变 , 射频 腔 体 的 品质 因数 相对 于 固定 频 
率 的 腔 体 要 低 得 多 ， 导 致 其 加 速 电压 相对 较 低 ; 男 一 方面 ， 超 导 同 步 回旋 加 速 
工作 在 ST 左右 的 强 磁 场 中 。 这 两 者 导致 超 导 同 步 回旋 加 速 器 的 轨道 圈 间 距 很 小 ， 
很 难 采 用 传统 的 静电 偏转 器 作为 引出 方式 ”; 同时 很 强 的 主 磁场 会 造成 负离子 的 
剥离 损失 , 因此 和 剥离 引出 也 不 适 和 应 用 于 超 导 同 步 回旋 加 速 器 ,， 超 导 同 步 回 旋 加 

器 一 般 都 只 能 加 速 正 离子 。 再 生 引 出 是 上 世纪 五 十 年 代 提 出 的 针对 同步 回旋 加 

速 器 的 一 种 引出 方案 ， 这 种 方案 最 终 被 Couteur Le 发 展 完善 ， 并 且 成 功 应 用 
于 IBA 230 MeV 质子 超 导 同 步 回 旋 加 速 器 和 Mevion S250i 质子 加 速 器 ”。 本 
文 将 针对 一 台 230 MeV 应 用 于 质子 治疗 的 超 导 同 步 回 旋 加 速 器 ， 对 其 再 生 引 出 
系统 进行 了 系统 的 动力 学 研究 ， 并 给 出 再 生 器 的 设计 。 
由 于 现 有 的 回旋 加 速 器 设计 程序 并 没有 针对 同步 回旋 加 速 器 再 生 引 出 系统 
的 设计 功能 ， 故 开发 了 能 够 进行 再 生 引 出 系统 设计 的 程序 。 同 时 引出 系统 的 设计 
和 主 磁场 紧密 相关 ， 因此 在 文章 的 第 一 部 分 , 我 们 将 介绍 这 台 超 导 同 步 回 旋 加 速 
器 的 主 磁铁 设计 ， 以 及 相关 的 轨道 、 动 力学 特性 和 加 速 器 主要 参数 。 文 章 第 二 部 
分 将 在 此 动力 学 程序 和 主 磁场 的 基础 上 , 提出 针对 这 台 加 速 器 的 再 生 引 出 系统 的 
设计 ,结合 轨道 动力 学 ,优化 再 生 器 的 结构 设计 ， 最 后 给 出 满足 要 求 的 引出 系统 
设计 方案 。 


1 主 磁 场 设计 
超 导 同 步 回旋 加 速 器 的 主 磁场 与 经 典 回旋 加 速 器 类 似 , 磁场 关于 中 心 旋 转 对 


称 性 ， 同 时 关于 中 心平 面 z=0 对 称 ， 磁 极 结构 比 等 时 性 回旋 加 速 器 的 扇形 磁铁 简 
单 。 为 了 减 小 加 速 器 尺寸 ， 参 考 IBA 设计 的 230 MeV 超 导 同 步 回 旋 加 速 器 所 采 
用 的 最 大 磁场 强度 ,该 同步 回旋 加 速 器 的 最 大 磁场 确定 为 5.71T, 位 于 磁铁 中 心 。 


根据 带电 粒子 在 旋转 对 和 
的 磁场 从 中 心 到 磁极 边缘 随 半径 应 逐 
考虑 到 磁极 的 利用 效率 ， 在 质子 主要 区 域 n<<1。 
速 电 压 一 般 较 低 ， 质 子 在 加 速 器 中 运动 的 圈 数 较 多 ， 所 以 受到 内 
在 主 磁场 设计 中 , 应 充分 考虑 共振 对 运动 稳定 
共振 线 。 本文 在 设计 主 磁场 时 考虑 了 五 阶 以 内 


过 不 可 避免 的 
响 。 当 质子 接近 
数 ， 以 最 小 化 却 


尔 磁 场 中 横 


向 运动 稳定 性 的 要 求 外 , 超 导 同 步 回 旋 加 速 器 
簿 渐 减 小 ,其 磁场 降落 指数 n 必须 满足 0<n<1。 


由 于 超 导 同 步 回 旋 加 速 器 的 加 
< 振 的 影响 较 大 。 


共振 线 时 调整 磁极 面 使 磁场 迅速 
振 对 质子 运动 的 影响 。 


示 。 主 磁铁 主要 参数 见 表 1。 


性 的 影响 号 使 质子 尽 可 能 


经 过 优化 ， 主 磁铁 1/4 截面 结构 如 图 


快 得 通 
共振 对 质子 运动 的 影 
过 共振 线 的 圈 
1 所 


束 减 小 , 减 小 质子 通 


表 1 主 磁铁 主要 参数 


Name 


Pole radius/mm 


Magnet yoke external radius/mm 


Magnet yoke internal radius/mm 


Height/mm 


Cover plate thickness/mm 


Magnetic pole gap/mm 


NI/A 


Weight/Ton 


Numerical value 


530 


1100 


750 


1360 


330 


128-68 


547x3000 


34.61 


图 1. 主 磁铁 1/4 截面 结构 图 ， 单 位 为 c 


图 2 为 通过 CST 模拟 ， 得 到 的 主 磁铁 中 心平 面 平均 磁场 及 磁场 降落 指数 n 
随 半 径 的 变化 曲线 图 。 磁场 沿 半径 的 变化 主要 通过 优化 磁极 面 的 结构 实现 。 在 磁 
铁 中 心 , 磁场 为 5.71T; 随 着 半径 增加 , 磁场 单调 减 小 , 在 束 流 引 出 半径 (r=44cm) 
处 ， 磁 场 降低 至 5.32T。 磁 场 降 落 指 数 随 半径 有 3 次 快速 增加 : 在 半径 8cm 内 ， 
磁场 降落 指数 快速 从 0 增加 到 0.017， 以 增加 质子 束 在 垂直 方向 获得 的 聚焦 力 ， 
从 而 提高 质子 束 的 俘获 效率 ; 在 半径 26cm 至 37cm 范围 之 内 , n 迅速 从 0.026 增 
加 至 0.065， 目 的 是 使 质子 迅速 通过 v=1/5 和 vz=1/4 共振 线 ; 磁极 边缘 处 磁场 降 
落 指数 迅速 增加 是 由 于 漏 磁 导 致 。 图 3 为 质子 在 超 导 同 步 回 旋 加 速 器 中 的 工作 路 
径 , 其 中 斜 线 表示 二 维 耦 合共 振 , 横 线 表示 一 维 共振 , 红 点 表示 从 半径 0 至 44cm 
范围 内 的 工作 点 变化 。 半 径 大 于 44cm 的 区 域 为 引出 区 ， 因 此 半径 大 于 44cm 的 
区 域 不 需要 考虑 。 考 虑 到 增加 再 生 器 后 产生 随 半径 增加 的 磁场 后 ， 会 使 w 增 大 、 
v 减 小 图 3 最 后 一 个 红 点 会 从 vs=1/3 共振 线 离开 并 移动 至 w=L/4 附近 ， 因 此 质 
子 主要 运动 区 域内 ，vz =1/4 共振 线 对 质子 的 运动 影响 很 小 。 由 图 3 可 知 质子 在 
加 速 器 中 受到 耦合 共振 的 影响 较 小 ， 但 会 受到 一 维 高 阶 共振 的 影响 。vz=L/3 共振 
在 磁极 边缘 区 域 ， 在 质子 主要 运动 区 域 之 外 ， 质 子 会 迅速 通过 ， 对 质子 运动 影响 
BUS; vz=1/4 和 vs=1/5 共振 可 能 对 质子 束 的 运动 产生 影响 ， 因 此 当 质 子 运 动 至 此 
共振 区 域 时 应 降低 磁场 使 质子 尽快 通过 ; 更 高 阶 的 共振 对 质子 束 的 影响 较 小 ,， 设 
计 主 磁铁 时 不 需要 考虑 。 当 然 ， 质 子 在 磁场 中 运动 的 详细 情况 ， 还 依赖 于 对 质子 
运动 的 详细 动力 学 分 析 。 
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图 2 中 心平 面 平均 磁场 和 磁场 降落 指数 n 随 半径 变化 图 


图 3 质子 工作 路 径 


2 横 、 纵 向 动力 学 分 析 
质子 在 电磁 场 中 会 受到 洛 伦 效力 ， 其 运动 满足 牛顿 第 二 定律 ， 即 : 
<P =qE+q(0xB) (1) 
其 中 , REP = ymoz 为 质子 的 机 械 动 量 ，? 为 质子 的 速度 , y = 1/V1 一 (v1c) 为 
相对 论 因子 ，c 为 真空 中 的 光速 ，mo 为 质子 静止 质量 。B 为 质子 所 在 位 置 的 电场 


强度 矢量 ，B 为 质子 所 在 位 置 的 磁感应 强度 矢量 ，g 表 示 质 子 电荷 量 。 


虽然 目前 的 三 维 电 磁场 计算 软件 可 以 得 到 精度 很 高 的 三 维 磁感应 强度 分 布 ， 
但 是 考虑 到 实际 磁场 测量 中 只 能 获得 磁极 中 心平 面 上 的 二 维 磁场 分 布 , 所 以 在 回 
旋 加 速 器 及 同步 回旋 加 速 器 的 轨道 计算 中 , 通常 都 是 根据 磁场 的 对 称 性 ， 利 用 中 
心平 面 的 二 维 磁场 分 布 ， 来 获得 磁场 在 三 维 空间 的 分 布 ”， 


Oo 10 1# 上 
十 十 B (r,0)+0 
ôr ror r OO” OES) (2) 


B_(r,0,z)= B.(r,0) La 


B(1,6,2)= z ČB, T, +0) 
or (3) 


By (1.8.2) == 5 B.(r.0)+ 02" r 

HPB, (r, 9) 为 中 心平 面 上 的 二 维 磁 场 分 布 。” 9,z 为 粒子 在 以 磁极 旋转 对 称 
轴 为 z 轴 ， 以 磁铁 中 心 为 原点 的 柱 坐 标 系 中 的 坐标 位 置 。 以 上 三 维 磁场 表达 式 只 
保留 了 偏离 中 心平 面 位 置 z 的 平方 项 ， 如 果 需 要 更 高 精度 ， 可 以 采用 类 似 的 方法 
展开 到 z 的 更 高 阶 。 本文 设计 中 ,考虑 到 在 引出 区 域 , 非 线 性 效应 的 影响 , 在 粒子 
轨道 计算 程序 中 ， 磁 场 展 开 至 z 的 平方 项 足够 使 用 。 

同步 回旋 加 速 器 中 , 粒子 只 在 加 速 间 际 受 到 电场 的 作用 。 该 电场 可 以 通过 谐 
振 腔 的 数值 计算 来 得 到 中 心平 面 上 的 电场 分 布 , 然后 利用 与 磁场 类 似 的 方法 , 得 
到 电场 的 三 维 分 布 。 


2.1 横向 动力 学 分 析 

在 不 考虑 引出 系统 影响 的 情况 下 ， 同 步 回 旋 加 速 器 中 粒子 的 轨道 比较 简单 ， 
其 静态 平衡 轨道 为 加 ， 非 参考 粒子 围绕 静态 平衡 轨道 振荡 , 振荡 频率 由 磁场 降落 
指数 决定 。 在 线性 近似 下 ， 其 轨道 满足 : 


2 
2 <+nz=0 (5) 
d0’ 
d’r 
pIE n)r =0 (6) 


其 中 z 表 示 粒 子 相 对 于 中 心平 面 的 位 移 ，r 表 示 粒 子 相 对 于 参考 粒子 的 位 移 ， 
0 表示 方位 角 ，n 为 磁场 降落 指数 。 通 过 方程 (5) (6) 可 知 非 参考 粒子 在 r、z 向 的 振 
荡 频 率 分 别 为 V1 一 n 和 Vn。 

周期 系统 的 周期 传输 矩阵 为 : 


JK (8) 
-JK sin(2x VK) cos(2x VK) 
p 


. cos(2z /KR ) 
| COS u ed 


—ysin 4 cos 


为 进行 多 粒子 模拟 分 析 , 横向 匹配 束 流 需 通过 周期 传输 矩阵 (8) 计 算 , EREC) 
表示 径 向 传输 矩阵 时 K=1-n, 表示 垂直 方向 传输 矩阵 时 K=n。 不 难得 出 B=1/y,，B = 
p/VK， 结 合 磁场 参数 ， 可 以 计算 半径 p 处 横向 匹配 柬 流 的 相 椭圆 参数 满足 : 


Br=tmax/T max=p/V1— N, B2=Zmax/Z max=p/Vne 

我 们 利用 自己 编写 的 程序 ， 对 粒子 运动 轨道 进行 计算 。 程序 采用 四 阶 龙 格 库 
塔 算法 ， 以 时 间 为 自 变 量 ， 磁 场 采用 式 (2)、(3)、(4) 展开 方式 获得 ，Bz(1,0) 关 
于 六 6 的 一 阶 和 二 阶 导数 通过 差分 方式 求 得 , 角 向 和 径 向 步 长 分 别 取 1 度 和 lcm， 
如 果 计 算 步 长 过 小 或 者 过 大 都 会 产生 较 大 误差 。 

粒子 自由 振荡 频率 随 半径 变化 如 图 4 所 示 , 图 中 点 线 表 示 的 振荡 频率 由 轨道 
计算 程序 计算 得 出 ,计算 时 先 确定 质子 静态 平衡 轨道 ,然后 在 静态 平衡 轨道 的 基 
础 上 改变 质子 横向 位 移 或 速度 方向 , 分 别 计算 质子 径 向 和 垂直 方向 的 振荡 频率 随 
半径 的 变化 关系 ; 实 线 表 示 的 振荡 频率 是 对 中 心平 面 平均 磁场 随 半 径 变 化 的 数据 
采用 差分 算法 直接 计算 得 出 。 由 图 可 知 轨道 积分 所 获得 的 横向 自由 振荡 频率 与 理 
论 值 符合 得 非常 好 , 说 明 程 序 计 算 精 度 满 足 要 求 ,虽然 磁场 展开 到 了 z 的 二 次 项 ， 
但 由 于 磁场 的 对 称 性 , 粒子 横向 速度 相对 角 向 速度 要 小 得 多 , 所 以 非 线 性 效应 并 
不 明显 。 
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图 4 质子 振荡 频率 随 半 径 变化 


2.2 纵向 动力 学 分 析 

纵向 动力 学 计算 时 , 我 们 假定 加 速 间隙 宽度 为 I cm， 加 速 间隙 之 间 的 电场 为 
均匀 分 布 ， 并 随时 间 按 下 式 变化 : 

E(t) = Eocos(wrs (Dt + ps)ëg (7) 
其 中 ，Bo 为 加 速 间 隙 的 电场 强度 峰值 ， 其 数值 上 等 于 加 速 电 压 峰 值 与 加 速 间 

隙 宽度 的 比值 。wrr 为 加 速 电 场 的 角 频 率 ，%s。 为 粒子 加 速 电场 初始 相位 ， 即 参考 
质子 平衡 相位 。 由 于 整个 加 速 过 程 中 相位 滑 移 系数 =ac-1/y* 恒 大 于 0, 则 平衡 相 
位 gs 必须 大 于 0. 

假定 参考 质子 同步 相位 为 30 E, 加速 电压 为 20kV, 使 参考 质子 经 过 加 速 
隙 获得 恒定 的 能 量 增益 。 编程 时 以 时 间 t 作 为 独立 变量 , 计算 的 时 间 步 长 在 107s 
至 10-10s 之 间 ， 步 长 为 非 恒定 值 。 通 过 计算 参考 质子 每 次 达到 加 速 间隙 中 心 的 时 
间 , 然后 参考 质子 从 第 一 次 达到 加 速 间隙 中 心 直 至 最 后 引出 时 加 速 电 场 的 相位 应 
依次 为 W6、pi/6+x、pi/6+2x……， 如 此 便 可 获得 电场 相位 随时 间 的 变化 关系 。 程 
序 计算 时 不 能 够 直接 获取 参考 质子 刚好 到 达 加 速 间 际 中 心 的 时 间 , 实际 计算 时 先 
获取 质子 经 过 加 速 间 隙 中心 的 前 一 个 时 间 点 和 后 一 个 时 间 点 , 然后 通过 线性 插值 
计算 参考 质子 刚好 到 达 加 速 间隙 中 心 的 时 间 。 通 过 式 (7) 可 知 电场 相位 随时 间 变 
化 的 一 阶 导 数 就 是 电场 频率 随时 间 的 变化 关系 ,参考 质子 运动 频率 和 经 过 加 速 
隙 时 的 电场 相位 随时 间 的 变化 如 图 5 所 示 , 图 5 显示 参考 质子 从 半径 5cm 加 速 到 
49.5cm 的 运动 过 程 ， 半 径 超 过 49.5cm 后 质子 基本 已 经 被 引出 因此 不 需要 考虑 后 
续 运 动 。 由 图 可 知 参考 质子 从 半径 5cm 加 速 到 49.5cm 需 要 约 0.106ms， 高 频频 率 
从 85.86MHz 减 小 至 60.61MHz， 电 场 相 位 从 pi/6 增加 至 pi6+15249r， 一 共 旋 转 了 
7624.5 图。 
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图 5. 参考 轨道 质子 运动 频率 随 半 径 和 时 间 变 化 关系 


一 一 初 相位 -29° 
-一 初 相位 65° 


初 相 A / 
| if 


Kinetic energy difference/MeV 


-60 -40 -20 0 20 40 60 80 
degree/° 


图 6, 相对 动能 随 相 移 的 变化 

图 6 为 5 个 质子 分 别 以 初 相位 30 度 、0 度 、-29 度 、65 度 和 -35 度 进 入 加 速 
电场 后 与 参考 质子 的 动能 差 值 随 相位 的 变化 关系 ， 由 图 可 知 以 参考 相位 30 度 进 
入 加 速 电 场 的 质子 动能 偏 移 量 几乎 为 0， 以 相位 0 度 进 入 加 速 电场 的 质子 虽然 偏 
离 平衡 相位 ， 但 能 随 参 考 质子 加 速 至 最 后 ， 整 个 加 速 过 程 中 其 相位 围绕 平衡 相位 
30 度 振荡 并 加 速 至 最 后 ， 最 大 动能 偏 移 量 大 约 为 0.7MeV; 以 相位 -29 度 进入 加 
速 电场 的 质子 其 相位 围绕 平衡 相位 30 度 振荡 并 随 参考 质子 加 速 至 最 后 ， 最 大 动 
能 偏 移 量 大 约 为 1.05MeV; 进入 加 速 电场 的 相位 分 别 为 65 度 和 -35 度 的 两 个 质 
子 ,最 终 会 进入 减速 区 域 并 损失 掉 , 说 明 处 于 相位 稳定 区 域 之 外 。 经 过 以 上 分 析 ， 
不 难看 出 非 参 考 质子 进入 加 速 电 场 的 相位 与 参考 相位 差 越 大 , 则 运动 过 程 中 非 参 
考 质 子 与 参考 质子 的 动能 离散 程度 越 大 。 由 图 可 知 当 平 衡 相位 为 30 度 时 相位 稳 
定 区 域 在 -30 度 至 60 度 之 间 ， 只 有 进入 相位 稳定 区 域 的 质子 才 会 随 参考 质子 加 


速 至 最 后 ， 没 有 进入 相位 稳定 区 域 的 质子 则 最 终 会 丢失 。 
多 粒子 模拟 的 相 空 间 随 加 速 圈 数 的 变化 情况 如 图 7 所 示 。 考虑 到 引出 时 要 求 

径 向 圈 间 距 大 于 1.3cm， 同 时 加 速 腔 对 束 团 垂直 方向 振幅 的 约束 ， 取 初始 束 团 径 
向 尺寸 为 1.4cm， 垂 直方 向 束 团 尺寸 为 2cm。 束 团 从 半径 5cm 处 开始 进入 加 速 间 
阶 ， 初 始 横 向 相 空 间 满足 : Br=rmax/r max0.0503，Bz=Zmax/z max~0.4226。 束 团 进入 
加 速 电场 的 相位 为 5 度 至 55 度 之 间 ， 初 始 束 团 的 动量 散 度 为 5% (a) 为 初始 时 
刻 半径 5cm 处 径 向 、 垂 直方 向 (横向 ) 的 匹配 相 空 间 分 布 以 及 纵向 相 空间 分 布 。(b) 
为 束 团 运动 至 22.9cm 运 动 2000 圈 ) 处 横 纵 向 相 空 间 分 布 。 (c) 为 束 团 运动 至 
33cm (运动 4000 Hel) 处 横 纵 向 相 空 间 分 布 。(d) 为 束 团 运动 至 引出 区 域 半 径 45cm 
(运动 6700 FA) 处 横 纵 向 相 空 间 分 布 。 由 图 可 知 质子 束 在 加 速 过 程 中 横向 相 空 
间 始 终 是 一 个 正 椭圆 , 这 是 由 于 束 团 在 同步 回旋 加 速 器 中 始终 受到 均匀 聚焦 力 的 
作用 ; 束 团 径 向 束 流 尺寸 变化 不 大 ， 这 是 由 于 整个 运动 过 程 中 径 向 振荡 频率 v, 变 
化 很 小 ， 垂直 方向 束 流 尺 寸 从 起 始 的 2cm 减 小 至 1.1cm， 这 是 由 于 整个 运动 过 程 
中 垂直 方向 的 振荡 频率 w 变 化 较 大 ; 整个 加 速 过 程 束 团 的 横向 夹 角 都 在 不 断 减 小 ， 
径 向 和 垂直 方向 的 最 大 夹 角 从 起 始 的 139mrad 和 23.66mrad 减 小 至 14.84mrad 和 
4.98mrad; 质子 束 纵向 动量 散 度 随 运动 圈 数 增加 而 减 小 ， 从 起 始 最 大 5%o 减 小 至 
1.8%o。， 整 个 运动 过 程 中 所 有 质子 的 相位 始终 围绕 平衡 相位 30 度 振荡 ， 纵 向 动量 
散 度 围绕 0 振荡 。 整 个 运动 过 程 中 没有 质子 丢失 。 
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图 7. 多 粒子 的 相 空间 随 运 动 圈 数 变化 


数 ， 见 表 2。 


基于 以 上 分 析 ， 最 终 可 以 确定 230 MeV 质子 超 导 同 步 回旋 加 速 器 的 主要 


Wp 


表 2 230 MeV 质子 超 导 同 步 回旋 加 速 器 主要 参数 


Name Numerical value 

Central magnetic field/T 5.71 

Extraction magnetic field/T 5.32 
Pole radius/mm 530 
Extracted beam energy/MeV 230 

RF frequency /MHz 85.8-60.5 

Acceleration voltage /kV 20 
Balanced phase/° 30 


3 再 生 引 出 系统 设计 
超 导 同步 回旋 加 速 器 加 速 电 场 幅 值 为 20kV， 轨 道 计算 程序 依据 本 文 的 设计 


参数 得 出 参考 质子 运动 圈 间 距 仅 为 40hm 左右 ,使 用 偏转 电极 作为 引出 器 ， 大 量 


质子 束 将 会 撞击 偏转 电极 造成 束 流 的 损失 ， 因 此 不 能 使 用 偏转 电极 作为 引出 器 。 
超 导 同 步 回 旋 加 速 器 一 般 使 用 再 生 器 引出 束 流 , 其 原理 是 在 磁极 边缘 垫 补 局 部 铁 
块 , 产生 沿 半径 非 线 性 增加 的 磁场 , 利用 v=1 共振 使 质子 经 过 再 生 器 后 轨道 中 心 
向 一 个 方向 移动 ， 同 时 经 向 圈 间 距 增 加 ， 垂 直方 向 仍然 维持 稳定 运动 。 一 般 认 为 
参考 质子 径 向 圈 间 距 大 于 1.3cm、 垂 直方 向 维持 稳定 运动 ， 质 子 加 速 全 230MeV 


时 具有 足够 的 径 癌 和 垂直 方向 稳定 区 域 ， 则 认为 再 生 器 设计 合理 。 


不 加 再 生 器 时 质子 在 加 速 过 程 中 会 遭遇 v=2v;z 共振 ， 受 到 w=2vz 共振 的 影响 


质子 径 向 运动 和 垂直 方向 运动 将 发 生 能 量 交 换 , 很 可 能 发 生 径 向 振荡 振幅 减 小 而 


垂直 方向 的 振荡 振幅 


兽 加 ， 从 而 造成 质子 束 的 损失 ,因此 再 生 器 的 放置 思路 应 是 


在 质子 进入 此 共振 之 前 将 质子 束 引出 加 速 嚣 外。 再生 器 的 设计 是 再 生 器 结构 调整 


与 动力 学 计算 之 间 反 复 迭 代 计 算 的 过 程 。 再 生 器 热 补 过 程 中 随 着 再 生 器 热 补 厚度 
的 增加 ， 再 生 器 产生 的 磁感应 强度 增加 ， 参 考 质子 的 运动 圈 间 距 也 随 之 增加 ; 当 


再 生 器 的 垫 补 厚度 超过 一 定 范围 以 后 , 参考 质子 的 运动 圈 间 距 虽 然 继续 增加 , 但 


其 垂直 方向 运动 振幅 也 会 迅速 增加 直至 撞击 磁极 而 损失 。 再 生 器 的 设计 应 当 满 足 
质子 符合 引出 要 求 的 同时 ,， 尽 可 能 增加 再 生 器 的 间 隐 使 尽 可 能 多 质子 通过 。 经 过 
不 断 和 迭代 计算 确定 再 生 器 最 终结 构 ， 图 8 为 再 生 器 最 终结 构 中 心 治 径 向 的 1/4 截 
面 图 和 再 生 器 半径 48.5cm 处 治 角 辐 截面 图 ， 红 色 轮 廓 表示 再 生 器 的 截面 图 ， 再 
生 器 结构 主要 参数 见 表 3。 
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图 8. 再 生 器 径 向 1/4 截面 图 和 再 生 器 角 向 截面 加 


表 3 再 生 器 主要 参数 
Radius/cm 36° 30° 24° 11° 
45-45.5 lcm 0.8cm 0.6cm 0 
45.5-49 1.05cm 1.4cm 1.35cm 0.06cm 
49-53 1.05cm Icm 1.1cm 0.06cm 
53-54.5 1.25cm 0.6cm 0.95cm 0 


结合 图 8 截面 图 和 表 3 再 生 器 主要 参数 可 知 再 生 器 沿 径 向 分 成 43cm-45.5cm、 
45.5cm-49cem.49cm-53cm 和 53cm-54.5cm 四 部 分 ,用 于 调整 磁场 随 半径 的 分 布 ; 
沿 角 向 可 分 成 11”、24”、30”、36” 由 部 分 , 用 于 减 小 再 生 器 中 心 区域 和 边缘 
区 域 磁 场 凹陷 引起 的 散 焦作 用 。 再 生 器 半径 为 45cm 至 45.5cm 范围 之 间 ， 从 远 
离 中 心平 面 至 靠近 中 心平 面 分 别 由 角 宽 36 度 lcm 厚 、 角 宽 30 FE 0.8cm 厚 和 角 
宽 24 EE 0.6cm 三 块 垫 补 铁 块 组 成 ， 再生 器 半径 为 45.5cm 至 49cm 范围 之 间 ， 
从 远离 中 心平 面 至 靠近 中 心平 面 分 别 由 角 宽 36 E 1.05cm 厚 、 角 宽 30 BE 1.4cm 
厚 、 角 宽 24 度 厚 1.1cm 和 和 角 宽 11 FE 0.06cm 厚 四 块 热 补 铁 块 组 成 ， 其 它 再 生 器 
模块 以 此 类 推 。 图 9 为 不 加 再 生 器 中 心平 面 平 均 磁场 随 半 径 分 布 、 单 独 再 生 器 产 
生 的 中 心 磁场 随 半径 分 布 以 及 磁极 增加 再 生 器 后 沿 再 生 器 的 中 心 磁 场 随 半径 变 


£. B| 10 为 半径 48.5cm 处 再 生 器 磁场 沿 方位 角 的 变化 。 
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图 9. 磁场 随 半径 分 布 
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图 10. 再 生 器 磁场 沿 角 向 分 布 
由 图 9 可 知 再 生 器 产生 了 一 个 0.8T 左右 的 峰值 磁场 , 同时 在 半径 为 42cm 和 
57cm 的 位 置 磁 场 为 -0.1T 左右 ， 总 磁场 的 峰值 约 为 5.98T。 由 图 10 可 知 再 生 器 在 
方位 角 158 度 和 202.6 度 的 位 置 磁场 为 -0.04T AA, 峰值 磁场 为 0.8T 左右 。 沿 径 
向 即将 进入 再 生 器 的 径 向 区 域内 和 刚 出 再 生 器 的 半径 区 域内 , 再 生 器 使 这 部 分 区 
域 的 磁场 减弱 ， 该 区 域 凹 陷 磁 场 对 质子 垂直 方向 的 运动 产生 散 焦作 用 ,因此 在 放 
置 再 生 器 时 需要 在 小 半径 区 域 垫 补 少量 铁 块 减 小 磁场 的 减弱 。 沿 角 向 再 生 器 的 边 
缘 和 中 心 区 域 也 会 产生 磁场 凹陷 , 该 部 分 凹陷 磁场 也 对 质子 垂直 方向 的 运动 也 产 
生 散 焦作 用 ,因此 再 生 器 的 角 向 边缘 需 呈 角 向 梯度 分 布 ， 同 时 在 中 心 区 域 热 补 一 


层 薄 铁 块 以 减弱 磁场 凹陷 。 
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图 11. 振荡 频率 随 参考 质子 动能 的 变化 

再 生 器 提供 随 半径 增加 的 磁场 ， 使 w 增 加 到 w=2/2 共振 的 同时 减 小 v,。 图 11 
为 主 磁铁 不 加 再 生 器 和 增加 再 生 器 时 w 和 2v, 随 参考 质子 动能 的 变化 关系 。 如 图 
所 示 不 增加 再 生 器 参考 质子 的 动能 在 248.2 MeV 时 达到 v,=2v,， 主 磁铁 增加 再 生 
器 后 v, 在 动能 为 237.1 MeV 时 增加 至 1，v; 为 0.278， 此 时 半径 为 43.7cm，237.1 
MeV 是 质子 可 以 达到 的 最 大 动能 。 

质子 束 的 引出 过 程 可 以 通过 图 12 所 示 径 向 相 点 图 理解 。 图 12 描述 的 是 再 生 
器 放置 在 方位 角 180 度 处 ， 质 子 运动 经 过 方位 角 180 REMY pE r 的 变化 关系 图 。 
12 中 描述 了 质子 动能 分 别 为 226.2 MeV、231.8MeV、234.0MeV、235.7 MeV 
和 237.0 MeV 时 径 向 稳定 区 域 的 边界 ， 以 及 质子 动能 为 231.8 MeV 时 两 个 不 稳 
定 轨道 。 图 12H r0 表示 质子 起 始 运动 半径 , 由 图 可 知 在 质子 动能 为 231.8 MeV 
时 径 向 稳定 区 域 为 3.5cm 左右 , 稳定 区 域 之 外 的 质子 径 向 运动 振幅 将 沿 着 渐 近 线 
呈 指 数 增加 。 此 外 观察 图 12 中 不 同 动能 的 质子 径 向 稳定 区 域 边界 可 知 随 着 质子 
动能 的 增加 ， 质 子 径 向 稳定 区 域 不 断 减 小 ， 当 ww 增加 至 1 时 质子 的 稳定 区 域 减 小 
至 0， 此 后 所 有 质子 将 处 于 不 稳定 状态 ， 径 向 圈 间 距 迅速 增加 ， 因 此 w=1l 是 质子 
可 以 达到 的 最 大 动能 , 同时 最 终 引 出 的 质子 束 存在 一 定 的 能 量 展 宽 。 通过 观察 图 
12 可 知 质子 动能 为 231.8 MeV 时 的 径 向 稳定 区 域 足够 俘获 绝 大 部 分 内 部 质子 束 ， 
所 以 该 再 生 器 可 以 加 速 质子 束 至 230 MeV. 
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图 12. 不 同 状态 的 7 个 粒子 径 向 相 图 
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图 13. 参考 质子 运行 最 后 6 圈 的 半径 随 方位 角 变化 

13 所 示 为 参考 质子 最 后 6 圈 径 向 轨迹 + 随 方位 角 6 的 变化 ,其 中 方位 角 的 
变化 范围 是 0 度 至 360 FE, 图 中 红色 矩形 框 表示 再 生 器 位 置 。 由 图 可 知 参考 质子 
在 方位 角 101 度 圈 间距 达到 1.5cm， 此 时 径 向 位 移 48.9cm， 磁 通道 便 可 以 在 此 后 
位 置 考虑 放置 。 如 果 选 择 使 用 无 源 磁 通道 , 无 源 磁 通道 将 被 主 磁场 磁化 并 对 其 周 
围 磁场 产生 一 定 的 影响 ,第 一 次 放置 磁 通道 后 还 需要 对 磁 通 道 的 位 置 进行 轻微 调 
整 ， 因 此 磁 通 道 的 放置 ， 也 是 位 置 调整 与 动力 学 计算 之 间 反 复 迭 代 计 算 的 过 程 。 
图 14 为 质子 从 半径 5cm 运动 至 满足 引出 条 件 时 垂直 方向 最 后 51 圈 的 轨迹 , 质子 
垂直 方向 初始 位 移 是 lem， 其 他 参数 都 与 参考 质子 参数 相同 。 整 个 运动 过 程 中 质 
子 在 半径 较 小 位 置 垂直 方向 运动 振幅 减 小 , 但 接近 再 生 器 后 垂直 方向 的 运动 振幅 
又 被 放大 ; 整个 运动 过 程 中 质子 垂直 方向 的 最 大 位 移 为 lcm， 质 子 垂直 方向 稳定 


z/cm 
i=] 
(=) 


0 5 10 15 20 25 30 35 40 45 50 55 
Tums 


图 14. 质子 最 后 51 圈 垂 直方 向 的 运动 随 运 动 圈 数 的 变化 
15 为 主 磁铁 增加 再 生 引 出 系统 后 多 粒子 模拟 的 相 空间 随 加 速 圈 数 的 变化 
情况 。 为 提高 程序 计算 速度 同时 不 影响 模拟 结果 ,选择 束 团 边 界 的 质子 进行 模拟 
计算 ; 半径 小 于 30cm 的 范围 ， 再 生 器 对 主 磁场 的 影响 很 小 可 以 忽略 不 计 ， 因 此 
选择 束 团 从 半径 30cm 开始 进入 加 速 电 场 。 束 团 的 径 向 和 垂直 方向 尺寸 与 主 磁铁 
多 粒子 模拟 一 样 取 1.4cm 和 2cm， 初 始 动量 散 度 为 5%。， 束 团 进入 加 速 电 场 时 的 
电场 相位 为 5 度 至 55 度 之 间 。 初始 横向 相 空间 参数 已 调整 至 匹配 。 图 中 (a) 为 束 
团 横向 和 纵向 初始 相 空 间 分 布 , (b) 为 束 团 运动 至 半径 34.2cm 处 (运动 1000 E) 
黄 纵 向 相 空 间 分 布 ，(c) 为 束 团 运动 至 半径 37.7cm 处 (运动 2000 圈 〉 横 纵向 相 
空间 分 布 ，(d) 为 束 团 运动 至 43cm 处 (运动 3000 圈 ) 动能 为 235.7MeV 时 横 纵 
向 相 空 间 分 布 ， 之 后 由 于 径 向 稳定 区 域 减 小 至 0， 束 团 圈 间 距 将 迅速 增加 。 由 图 
可 知 在 加 速 过程 中 相 空间 逐渐 畸变 , 这 是 由 于 主 磁铁 增加 再 生 器 后 磁场 对 束 团 的 
聚焦 效果 为 非 均匀 育 焦 ; 束 团 运 动 半径 增加 ， 相 空间 畸变 越 严 重 ,， 束 团 运动 半径 
大 于 40cm 后 相 空 间 畸 变 特别 明显 ; 束 团 从 半径 30cm 运动 至 43cm 的 过 程 中 ， 
纳 向 尺寸 和 夹 角 并 没有 发 生 显著 增加 ,纵向 动量 散 度 从 起 始 最 大 5%o 减 小 至 2.8%o， 
纵向 相位 宽度 所 增加 但 都 在 稳定 区 域 之 内 ; 整个 运动 过 程 中 质子 束 运 动 稳定 , H 
没有 质子 丢失 ， 说 明 束 团 可 以 稳定 加 速 至 235.7MeV。 结 合 以 上 分 析 ， 主 磁铁 增 


235.7MeV 。 


加 再 生 引出 系统 后 具有 足够 的 径 向 和 垂直 方向 稳定 区 域 以 加 速 内 部 质子 束 至 


图 15. 多 粒子 模拟 的 相 空间 随 运 动 圈 数 变化 情况 


16 为 初始 束 团 横向 相 空 间 分 布 ， 束 团 中 心 为 半径 43cm， 质 子 动能 皆 为 
235.7 MeV, 工 同 尺寸 分 别 为 1.4cm 和 1.0cm, z 向 的 尺寸 分 别 为 6mm 和 lcm， 初 
始 束 团 相 空 间 参 数 已 调整 至 匹配 。 假 定 磁 通 道 入 口 位 置 为 方位 角 105 FE, AEE 


48cm 至 50cm， 此 处 参考 质子 的 圈 间 距 为 1.7cm， 图 16 所 示 初 始 束 团 运动 至 磁 
通道 入 口 位 置 时 的 横向 相 空间 分 布 如 图 17 所 示 。 图 17 中 黑色 和 矩形 框 表示 磁 通 道 


入 口 的 位 置 ，r=0 表示 半径 49cm， 图 17 中 显示 的 是 初始 束 


romax=7mm 和 Zomax=3mm、 romax=7mm 和 Zomax=Smm、 romax=Smm 


团 相 空 间 分 别 以 
和 Zomax=3mm 为 


初始 参数 ， 最 后 计算 所 得 结果 。 由 图 可 知 束 团 在 再 生 器 的 影响 下 r 向 和 z 向 的 尺 
寸 和 运动 夹 角 均 被 放大 ; 束 团 rt 向 尺寸 都 基本 差不多 大 ,人 磁 通 道 入 口 位 置 处 r 向 


最 大 位 移 分 别 为 9.5mm、10.0mm 和 9.7mm; RA z 向 尺寸 放大 较为 明显 ，z 向 
最 大 振幅 分 别 为 9.5mm、13.9mm 和 8.0mm， 因 此 要 想 提 高 束 团 的 引出 效率 ， 约 
束 束 团 垂直 方向 的 尺寸 效果 较 明显 。 考 虑 到 再 生 器 的 间 隐 约 为 3.2cm， 束 团 运 动 


至 43cm 处 时 ，z 向 尺寸 最 好 应 维持 在 Smm 以 内 。 
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图 16， 初 始 束 团 横向 相 空 间 分 布 
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图 17， 束 团 运 动 至 磁 通道 入 口 处 横向 相 空间 分 布 


本 文 针 对 230MeV 质子 超 导 同 步 回 旋 加 速 器 开展 了 再 生 引 出 系统 的 设计 研 
究 ， 由 于 再 生 引 出 系统 与 主 磁铁 紧密 相关 ， 也 给 出 了 主 磁铁 的 设计 研究 。 现 有 的 
回旋 加 速 器 设计 程序 并 没有 针对 同步 回旋 加 速 器 再 生 引 出 系统 的 设计 功能 , 所 以 
本 文 开发 了 能 够 进行 同步 回旋 加 速 器 再 生 引 出 系统 设计 的 程序 。 本 文 设计 的 主 磁 
铁 极 面 结构 相对 简单 同时 关于 中 心平 面 旋转 对 称 , 目的 是 为 了 降低 工程 上 的 实现 
难度 。 动力 学 计算 结果 显示 中 心平 面 磁场 满足 径 向 和 垂直 方向 聚焦 条 件 ， 质 子 在 
加 速 过 程 中 迅速 通过 共振 区 域 可 以 正常 加 速 到 引出 区 域 表 明 主 磁铁 的 结构 设计 
满足 要 求 。 利 用 程序 跟踪 参考 质子 最 后 儿 圈 的 运动 状态 , 质子 最 后 圈 间 距 可 以 大 
于 1.5cm 同时 垂直 方向 的 运动 仍然 维持 稳定 。 本 文 设 计 的 主 磁铁 和 再 生 引 出 系统 
具有 足够 的 径 向 和 垂直 方向 稳定 区 域 ， 以 加 速 内 部 质子 束 至 动能 235.7 MeV; 但 
随 着 质子 束 进一步 向 外 运动 ,垂直 方向 的 运动 振幅 放大 较 大 , 根据 本 文 动力 学 计 
算 结 果 显示 ， 质 子 束 加 速 至 235.7 MeV 时 垂直 方向 的 运动 振幅 最 好 维持 在 5mm 


以 内 ; 动力 学 计算 结果 显示 磁 通 道 入 口 处 不 可 避 倪 会 造成 质子 束 损失 。 根据 再 生 


器 处 的 动力 学 模拟 结果 可 知 ， 通 过 限制 加 速 器 中 心 区 域 质 子 束 垂 直方 


可 的 振幅 ， 


可 以 减 小 质子 束 在 再 生 器 区 域 的 损失 , 进而 提高 引出 效率 , 同时 减 小 质子 束 碰撞 
加 速 器 结构 导致 的 活化 。 
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Design Study of 230 MeV Proton Superconducting 


Synchrocyclotron Regenerator extraction system” 
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Institute ofNuclear Science and Technology, Sichuan University, Chengdu 610064, 


China) 


Abstract: Superconducting synchrocyclotron has an ultra small structure, lower 
cost and is more suitable for use in proton therapy systems used in hospital 
environments. Compared with the current proton therapy system, the proton therapy 
system using the superconducting synchrocyclotron has a lower cost for single 
tumor treatment and has a huge application prospect, so it has been highly valued 
by many accelerator research and development institutions and enterprises 

Regeneration extraction system is the key system of superconducting 
synchrocyclotron, and it is also the key and difficult point in the design of 
superconducting synchrocyclotron. In this paper, a regeneration extraction 
system was designed for a 230MeV proton superconducting synchrocyclotron for 
proton therapy. Because the regeneration extraction system is closely related to 
the main magnetic field, the design of the main magnet is carried out in the 
first part of this paper, and the related dynamic calculation results are given. 
The existing cyclotron design program does not have the design function of 
synchrocyclotron regenerative extraction system, so this paper develops a program 


that can design the regenerative extraction system. The design results show that 


the design of 230 MeV proton superconducting synchrocyclotron regeneration 
extraction system meets the extraction requirements, which can provide a 
reference for the design and engineering construction of other systems of 
superconducting synchrocyclotron. 

Key words: superconducting synchrocyclotron; subminiature structure; design 


regenerator extraction system; design main magnet; 
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